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Abstract 

We give an algorithm for finding network encoding and decoding equations for error- 
free multicasting networks with multiple sources and sinks. The algorithm given is 
efficient (polynomial complexity) and works on any kind of network (acyclic, link 
cyclic, flow cyclic, or even in the presence of knots). The key idea will be the 
appropriate use of the delay (both natural and additional) during the encoding. The 
resulting code will always work with flnite delay with binary encoding coefficients. 

Keywords: Network codes, multicasting, cyclic networks, pipelining, delay. 

1 Introduction 

We consider centralized algorithms for designing the network coding equations for mul- 
ticast [T] in a network with directed edges. Early work in this area focused on acyclic 
networks (see |3], |1] and the references therein). Cyclic networks were considered in 
and in [6], and [H]. In [5] we developed, based on the Linear Information Flow (LIF) al- 
gorithm in |3] , the LIFE and LIFE-CYCLE algorithms. These will find a linear encoding 
(if one exists) in networks that are, respectively, link cyclic and simple flow cyclic (see [5] 
and Section [2] for definitions of different types of cyclicity) . 

Thus the state of the art prior to this paper is that centralized encoding can calculate, 
in polynomial time, a valid encoding for acyclic and simple cyclic networks. However, 
many networks occurring in practice may, in fact, contain complex cyclic structures. 

There are two contributions of this paper: 

1. We extend the algorithms in by applying Mason's formula tll],[T5] for signal 
propagation in a cyclic graph. The resulting new algorithm (LIFE*) will calculate a 
valid encoding for arbitrary cyclic and acyclic networks. The complexity is similar 
to that of the LIFE algorithm. 

2. We propose a simple binary encoding scheme that exploits the natural delay inherent 
in the network. This scheme is related to those proposed in [8], [11], [12], and [13], 
and also to the randomized version of the LIF algorithm [1]. 



1 



In Section |2] we give an overview of the necessary notation and previous results. Fur- 
ther, we will recall the two different notions of cyclicity and the two different types of flow 
cycles. Next, in Section |3] we present the LIFE* algorithm. Although the algorithm works 
over any field, we propose that the delay based binary encoding scheme demonstrated in 
Section |3] is particularly well suited for the LIFE* algorithm. Finally in Section |4] we 
discuss complexity issues and practical aspects. We have included an appendix, which 
makes up the bulk of the paper, with detailed examples that show how the algorithms 
work on the different kinds of networks. 

2 Notation and previous results 

The notations used in this paper follow and expand those used in 

Consider the network G, where G = (V, E) is a directed multigraph. V represents the 
set of nodes and E the set of edges or links, each of unit capacity. A pair of nodes may be 
connected by one or multiple edges. Let 5* = {si, 32, ... , Sh} C \^ be the set of unit rate 
information sources and let T = {^1,^2, ■ ■ ■ ,tr} C V he the set of sinks. We will assume 
that each source (synchronously) generates one symbol at each (discrete) time instant. 
Let X G Z denote the time. We will denote by crj(x) G F2 the symbol produced by source 
Si at time x. In case the given network has one unique source s that sends information at 
rate h we will simply create virtual sources Si, . . . ,Sh and link them to the actual source 
s. Note that most real life networks can be precisely or approximately described by a 
network G as outlined here. 

For each edge e G -E we will denote start{e) and end{e) the nodes at which e starts 
and ends, respectively. As usual, a path from node u to node v of length I is a sequence 
{ci E E : i = 1, . . . ,1} such that u = start{ei), v = end{ei), and start{ei+i) = end{ei), for 
i = 

The complete set of symbols ai{x), i = 1 . . . h will be called generation x. The object 
of the algorithm is to find an assignment of equations in such a way that each sink t G T 
at each time x + dt can complete the decoding of the whole generation x, where dt is a 
(finite) constant for t denoting the total delay associated to that sink. It is assumed that 
cri{x) = 0,i = 1 . . . h for any x < 0. 

We assume that the transmission of a symbol on each link e has a unit delay associated 
with it, that is to say, if a symbol s is being carried by edge e at time x, that symbol, 
once processed at node end{e), can be carried by an edge e' at time x + 1, for any edge 
e' with start{e') = end{e). Apart from this intrinsic delay, as we will see, the encoding 
process might assign extra delays in order to satisfy certain required conditions. 

Let D denote the linear delay operator. The way in which the operator works is as 
follows: 

D{a{x)) = Da{x) = a{x - 1) 

and it is extended by linearity. 

D{a{x) + ^(x)) = Da{x) + D^{x) = a{x - 1) + ^(x 

DV(x) = D{D'-^a{x)) = a{x - i) Vi G Z+ 

{D' + D^)a{x) = D'a{x) + DV(x) = a{x - i) + a{x 

Also, abusing notation, we extend the operator in order to work 
nents in the following way. 



-1) 
-j) 

with negative expo- 
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L)-V(x) = (t(x + 1) 

A flow path is simply a path from a source s to a sink t. We will assume that for 
each sink t & T there exists an edge disjoint set of flow paths /* = {/^*'*, i = 1, . . . h}. 
Following the notation of [3J we will call /* the flow for sink t. The minimal subgraph of 
G that contains all flows (and their associated nodes) will be called the flow path graph, 
and can be determined from G by a suitable polynomial algorithm. /^From the perspective 
of encoding as discussed in this paper, we ignore the issue of determining the flow path 
graph and the flows, and assume that the network is a flow path graph, that the flows are 
known, and that every edge is on some flow path. 

The results in pQ guarantee that in such a network, all the sinks can receive all the h 
input symbols produced by the h sources, and the results in j2] state that it can be done 
using linear coding on the network. In this paper we deal with linear coding, so each edge 
will be encoded using a linear encoding equation. 

Definition 1 A link cyclic network is a network where there exists a cyclic subset of 
edges, i. e., a set {ei, 62, . . . , e^, Cfc+i = ei} C E for some positive integer k such that 
end{ei) = start(ej+i) for 1 < i < k. The set of edges {ei, 62, ... , e^, e^+i = ei} is a link 
cycle. If no such cycle exists, the network is link acyclic or simply acyclic. 

Suppose e is an edge that lies on the flow path /*"*^. We will denote by ft^{e) the 
predecessor of edge e in that path. There is no ambiguity in the notation: e can lie on 
several flow paths for different tj, but not for different Si and the same tj, since all 
the flow paths arriving in tj are edge disjoint. Thus, once tj is fixed, e can only have one 
predecessor in the flow to tj, that is, ft^{e) is the predecessor of e in the only flow path 
arriving in tj that contains e. In the same way we denote by /-^(e) the successor of e in 
that flow path. 

Let T(e) C T denote the set of sinks t that use e in some flow path /*, and let 
-P(^) = {/i-(^) I ^ ^ ^(^)} denote the set of all predecessors of edge e. Note that all edges 
will have some predecessor (P(e) 7^ 0), except those with start{e) = Sj for i G {1, . . . , h}. 

We will introduce some extra notation to denote the temporal order induced in the 
edges by the flow paths in /. When two edges Ci and 62 lie on the same flow path and 
ei is the predecessor of 62 in that path, we will write ei -< 62. We will use transitivity 
to define relationships among other pairs of edges that lie on the same path but are not 
consecutive. We observe that in each path the relation -< defines a total order in the 
edges that form that path, since a path that contains cycles, that is, edges that satisfy 
Ci -< €2 -<■■■-< Cn -< Ci can be simplified by avoiding taking the trip around that cycle. 

Definition 2 A flow acyclic network is a network where the relation -< defines a partial 
ordering in E. If the relation does not define a partial ordering, the network is flow cyclic. 

Note that a flow cyclic network is always link cyclic, but the converse is not true. 

Definition 3 A simple flow cycle is a link cycle {ci, 62, ... , e^, e^+i = ei} such that for 
each i = 1, . . . ,k there exists a flow path /'^*'*' that implies Ci -< Cj+i. We observe here 
that must be at least two distinct flow paths traversing the cycle. 

Definition 4 A flow knot or simply a knot is formed by two or more simple flow cycles 
that share one or more edges. 
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Figure 1: A flow knot. 



Figure 1 shows an example of a knot. This particular knot forms part of the network 
in Example 4 in Figure 2. 

We illustrate these concepts with some examples, shown in Figure 2, in which we show 
the flow path graphs for some networks. To make it easier to see, each sink is represented 
with a different color and all the flow paths arriving in that sink are drawn in the same 
color. One can easily check that the first example is an acyclic network (it is one of the 
so-called combination networks), the second is link cyclic but flow acyclic while Examples 
3 and 4 are flow cyclic. Example 3 is the same cyclic network presented in |T], while 
Example 4 has been created by us to illustrate a network with a flow knot. 



In the whole process of encoding we will observe two basic principles. 

• Pass on information (PI) principle: [5] When encoding each edge e G -E, all 
the symbols carried by the incoming edges f^{e),t G T(e) will contribute. This 
is because assigning coefficient to any of those predecessors means that the flow 
carried on that predecessor is stopped at that point, meaning in turn that the actual 
flow path graph used is different from the one established initially, j^] 

• Old symbol removal (OSR) principle: [5] To avoid symbols circulating endlessly 
on each flow cycle, they should be removed at the entrance point by the node through 
which they entered the cycle. We will explain this in more detail in |3.2 

In what follows we will consider two ways of expressing the symbol Ve{x) G F2 to be 
transmitted on edge e at time x. 

The local encoding equation specifies the action of each node. The general local en- 
coding equation is 



where 7r{p, e) is a polynomial in the operator D which denotes the encoding, i. e. they 
express the additional delays introduced for each flow path at the node start{e). The 
term r(p, e) denotes a function of the operator D, which we will call the transfer function 
from ptoe and which accounts for the natural delay that, as previously noted, is inherent 
to the transmission on each edge. The way to compute this function will be explained in 

* Observe that this principle makes sense if and only if the flow path graph is "reasonably efficient" . It 
may require some effort to calculate a reasonably efficient flow path graph, but if that work has been done 
it is wise, in terms of complexity, to rely on the provided flow path graph. The connection between flow 
path graph calculation and network coding can offer interesting complexity trade-offs, but this discussion 
is beyond the scope of this paper. 



3 The LIFE* algorithm 




(1) 



peP(e) 
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Figure 2: Examples of different types of networks. 
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detail for each case. The process of determining the network code consists of finding a set 
{Ti{p,e),ee E,pe Pie)}. 

Remark: To simphfy the encoding, we propose to use only monomial 7r(p, e)'s, that 
is, vr(p, e) = for some ie{p) G Zq . Thus the encoding consists of simply adding the 

encoding vectors of all predecessor edges, each one artificially delayed for zero or more 
time units as necessary. The algorithm can be straightforwardly modified to allow all 
polynomials in ¥2{D) to be coefficients of the encoding combination. We have chosen 
the proposed scheme because it is simpler and because in this way we stick to the flow 
path graph that has been computed beforehand. Since there exists an encoding for any 
flow path graph, once a flow path graph has been computed we want to follow it and 
proflt from that in order to reduce the computational complexity of finding the encoding 
equations for each edge. With this simplification, the encoding takes the form 

v,ix)= J2 D^'^'^r{p,e)v,{x), (2) 

peP(e) 

where ieij)) ^ exponents that express the additional delays introduced at the 

node start (e). 

The exponents ie{p) do no depend on the time variable x. We remark that employing 
a time invariant encoding, as in this case, ensures that streams of symbols from each 
source can be pipelined through the network. 

The global encoding equation expresses the symbol Ve{x) to be transmitted on edge e 
in terms of the source symbols: 

h 

Ve{x) = a{e,i,y)ai{x -y), (3) 

where a{e, i, y) G F2 is the coefficient that specifies the infiuence of each source symbol 
on Ve{x). 

Since cTj(z) = for all z < 0, the sum over ?/ in ([s]) is finite, since the terms with y > x 
will all be 0. 

Another way of expressing the global encoding equation, more useful for our purposes, 
is by means of the field of rational functions ¥2(0), where D is the delay operator. 

h 

Ve{x)=J2FeAD)ai{x), (4) 
i=l 

where Ff.i{D) G ¥2{D). Again, all the monomials in with y > x will give null terms 
in the sum. 

The goal of the LIFE* algorithm is to determine the local encoding equations of each 
edge e in each flow path, and the global encoding equation of the last edge in each flow 
path. In fact, during the course of the algorithm, all the global encoding equations will 
be determined. 

In particular, the global encoding should be such that each sink t can, after a certain 
constant delay df, extract all the h information symbols belonging to a generation from 
the h sources. Each sink t receives h symbols, namely {vei{x) \ei G f^^'^,t = end{ei),i = 
1, . . . , h}. If the corresponding h global encoding equations are linearly independent over 
the fleld ¥2(0), the sink t can reconstruct the h input symbols of generation x at time 
X + df, by solving (in a simple way) the corresponding set of equations. 
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In the case of a flow acyclic network, all the edges in E can be visited according to 
the partial order -< induced by the flow paths, starting by the edges with no predecessor 
(e e -E I start{e) = Si for some i), and proceeding in such a way that an edge e will not be 
visited until all its predecessors (all p G -P(e)) have already been visited and their global 
encoding equations computed. If the global encoding equations of all predecessors of e 
are substituted into ([2]), we get the global encoding equation for e. 

When the network is flow cyclic the algorithm comes to a point at which no more edges 
can be visited in topological order, that is, the algorithm hits the flow cycle. Once the 
edges taking part in that cycle have been identified, the encoding of the whole cycle has 
to be treated as a whole. Here the old symbol removal principle is useful. An appropriate 
version of Mason's formula (see [H]) for the transfer function on a circuit will be the main 
tool to deal with it. We will explain this in detail in the second part of this section. 

Similar to the LIF and the LIFE algorithms, the new LIFE* algorithm proceeds by 
maintaining, through the iterations of the main loop, the following sets for each t E T: 

• A set Et G E , \Et\ = h, such that Et contains the most recently visited edge on each 
flow path in /*. 

• An h X h matrix Mj, in which the element i,j will be the coefficient of (Ji{x) in the 
global encoding equation of the j-th element of the subset Ef. 



\ Fe,„,(Z}) ■■■ Fe,„,(D) ) 

where Ct^i, . . . , Ct^h are the edges in Et. 

Through each step of the algorithm we will impose the full rank condition: The 

matrix Mt must have rank h, for all t E T. This condition is sufficient, but not necessary, 
for obtaining a valid network code. At the final step, for each t e T the set Et will be 
that of the h edges which arrive in t. If these edges carry the symbols rt^i{x), . . . ,rt^h{x) 
at time x, we have 

{ai{x),...,ah{x))Mt = {rt^i{x), . . . ,rt^h{x)) ■ 

By the full rank condition the matrix Mt is invertible, and the symbols {ai{x), . . . ,ah{x)) 
can be recovered from the received ones as 

{ai{x), . . . ,ah{x)) = {rt^i{x),. . . ,rt,h{x)) Mt~'^. 

We will now explain how the algorithm will proceed. 



3.1 The flow acyclic parts of the network 

Whenever there is an edge e that can be visited in topological order (that is to say, one 
for which all the predecessors have been already visited), the algorithm will proceed to 
update the sets Et and Vt in the following manner: 

Et := {Et \ {/1(e)}} U {e}, for each t e T(e) 
We consider the general local encoding equation for edge e 
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where ie{p) G Zq are the unknowns and represent the extra delay added at edge e to 
maintain the full rank invariant. 

In the acyclic case with unit link delay, t(p, e) = D, G -P(e). If the natural link 
delay is any other thing, even not necessarily the same for all links, it will be straightfor- 
ward to adapt the corresponding equations. The local encoding formula ^ in the acyclic 
case with unit link delay takes the form 

v^{x) = D'^^^^Dvpix) = D'^^^^+hp{x) (5) 

When the global encoding equations for all the Vp{x), already determined in the pre- 
vious steps of the algorithm, are substituted in the above expression, we will have the 
global encoding equation for e. 

Further, for each t G T(e), replace in each matrix Mj the column corresponding to the 
encoding equation of /^(e) with the new one Ve{x). 

The unknowns will be chosen to have values that satisfy the full rank invariant for 
all the matrices Mt,t G T(e), that have been updated following the update of the corre- 
sponding set Et. 



Conjecture 1 There exists a finite value I, that depends on the graph and in particular 
on the set T{e), so that some set {ieip) '■ P £ -P(e) (^iT-d ie{p) < I}, when applied to 
and pi), will satisfy the full rank condition. 



This conjecture is in a way similar to Lemma 6 in [3] . There the field size required to 
guarantee the full rank invariant is proven to be |T(e)|, that is to say, each coefficient of 
the linear combination can be chosen among |T(e)| possibilities. If we assume / = |T(e)|, 
the coefficients of the combinations in (12) can be chosen also among |T(e)| different 
possibilities, namely = 1, D, D^, . . . , £)l^(^)l~i. Nevertheless, the proof cannot use the 
Linear Algebra arguments used there because the set with which we are working is not a 
vector space. 

The conjecture is further supported by Theorems 1 and 3 in [13], and by software 
simulations for random networks that we have carried out. We omit the details of these 
simulations. 

For networks with random structure most encodings need no extra delay at all, and 
in the few remaining cases a delay of one unit in one of the incoming paths is enough 
to solve the problem for most of them. This observation also tells us that in order 
to find the encoding for each edge, which means finding the delay exponents ie{p), an 
efficient approach will be to use a greedy algorithm (see [3]) that starts by considering 
^e(p) = 0, Vp G -P(e), checks if the full rank condition is satisfied, and in case it is not, 
proceeds to increment the delay exponents one by one until a solution is found. According 
to our simulations, the average number of tries needed to find the solution for each edge 
will be very low. 
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3.2 Dealing with flow cycles 

When the algorithm encounters a flow cycle the set of edges forming part of the cycle has 
to be computed. Let us call Ce the set of edges and Cy the set of nodes that are ends of 
those edges in the cycle. 

Let us call P{C) the set of predecessor edges of the cycle, 

P(C) = {ee E\CE\ee P(e') for some e' G Ce}- 

We assume that the encoding equations for all the edges e e P{C) have been already 
determined in previous steps of the algorithm. The goal of this step of the algorithm will 
be finding at once the encoding equation of all the edges in C^;. In a sense it is as if the 
cycle as a whole is being treated as a kind of 'superedge'. All the individual edges in 
it will have basically the same structure of equation, which will be a combination of the 
equations of the edges in P{C), that is 

vc{x)= J2 D'^^''^r{p,C)v^{x) (6) 
peP{C) 

Here r(p, C) is a notation with which we simply mean a transfer function that will 
have to be computed separately for each particular edge in Ce- Thus, each edge in Ce 
will have a slightly different version of that basic structure due to the fact that they lay 
in different parts of the cycle and will consequently observe the incoming equations with 
different delay. We will distinguish two cases, namely, when the flow cycle is simple, or 
when it is a knot. The distinction will be made just for the clarity of explanation, since 
the simple case is just a particular case of the knot case. 

Once the exact encoding equation has been determined for all the edges in C^, the full 
rank invariant has to be checked only for the last edges in the cycle for each flow path, 
that is to say, for each t with (1 Ce ^ ^ the corresponding full rank condition must be 
satisfied by the edge e e HCe such that /l^(e) ^ Ce- 

3.2.1 The simple flow cycle case 

Suppose Ce = {ei, 62, ... , e^} with end{ei) = start{ei-i) for i = 1, . . . , k—1 and end{ek) = 
start{ei). 

Let us consider the local encoding equation of the cycle ([6]). We show now how to use 
the Old symbol removal principle. 

We will focus on a certain edge in Ce, for instance Ci. Suppose the local (and global) 
encoding equations of Cfc, the predecessor of ei in the cycle, have been determined exactly. 

p£P{C) 

If ei had no predecessor outside the cycle, that is to say P(ei) = {e^}, then the en- 
coding equation of ei would simply be Vej^{x) = Dve^{x) = J2peP{C) ek)Dvp{x), 
which means r(p, Ci) = Dt{p, e^) G P{C). In the same way it is clear in general that 
if p e P(C) \P(ei), then r(p,ei) = Dr^p^ek)- 

Now suppose that there is one particular edge pi G P{C) which is the unique prede- 
cessor of ei not in Ce-, this implies r(pi,ei) = D. Not removing at that point the old 
contribution of that predecessor would mean that the new contribution would mix with 
the old ones on each loop of the cycle and would keep circulating forever. In order to 
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avoid this we want to remove the old contribution that came on edge pi k time instants 
ago (where k is obviously the length of the cycle) and contribute to the circulation in the 
cycle with only the newest symbol coming on pi. This is done as follows 

Ve,{x) = D ■ Ve.ix) + D'^^P^W{pi, d) [-Vp,{x ~ k) + Vp,{x)] 
= D-Ve^ (x) + D'C'-P^W [-Vp^ {x-k)+ Vp, (x)] 

Here we have used the minus operator (— ), despite all the operations are always on the 
binary field, in order to emphasize which symbols are being removed from the circulation. 

In general, if e G C^; has several predecessors not lying in the cycle, the local encoding 
equation of e in terms of the predecessors of e takes the form 

Ve{x) = D ■ Vp^(e){x) + EpeP(e)nP(C) D'cip)T{j>, e) [-Vp{x -k)+ Vp{x)] 
= D ■ Vp^(e){x) + EpeP(e)nP(c) D'c^^~^D [-Vp{x - k) + Vp{x)] 

where pc{,e) = P(e) fl Ce- 

The result of doing this at the entrance in the cycle of each predecessor p G P{C) is 
that only one 'instance' of the symbols carried by each p G P{C) will be circulating on 
each edge e E Ce- It is easy to see that the transfer function from each predecessor of the 
cycle to each edge in the cycle will be r(p, e) = d{p, e) Vp G P{C), e G Ce, where d{p, e) 
is the 'distance' measured in number of edges in the cycle that lay between end{p) and 
end{e). 

To summarize, the local encoding equation of each edge e G C^; in terms of the 
predecessors of the cycle is 

v^{x) = J2 D'^^v^+'^^P^^^Vp{x) 
peP{c) 

where again, the only unknowns are the values ic{p) £ '^o- 

An example of the use of this procedure can be found in the Appendix when encoding 
the network of Example 3. 

3.2.2 The knot case 

Suppose the set C^; is not just a simple flow cycle but forms a knot. 

Again the local encoding equation for the whole knot will share a common structure 

vcix)= D'^^''^r{p,C)vp{x) 

P<^P(C) 

Once more the idea is to treat the whole knot as a kind of 'superedge'. 

Here the main idea is the same as before: the old symbols must be removed from the 
circulation. In order to do it one needs to know how those arrive at each edge of the knot, 
and for this we need as a tool Mason's formula (see [Il],[l5]) for the computation of the 
transfer function on a cyclic circuit. 

We apply Mason's formula to the directed line graph associated with C^; in the fol- 
lowing way: Two edges e' and e in C^; are considered adjacent (and an arc starting in e' 
and ending in e will be drawn in the line graph) if and only if there exists a flow /* for 
some t G T such that /i_(e) = e'. 

For each symbol entering the knot we have to compute the corresponding transfer 
function over all the edges in C^;. For this we will consider the entrance point of the 
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symbol (the edge at which that symbol enters) and the exit point (that is, the edge at 
which we want the transfer function of that symbol), and will apply Mason's formula 
between these points in the line graph above mentioned and using the delay operator D 
as the branch gain (see |14|) of each edge. Again, each branch gain can be taken to be 
whatever function of D models best the actual behavior of the transmission on that edge 
and the corresponding equations can be adapted consequently. 

We show in detail how to compute the transfer functions by means of Example 4 in 
the Appendix. 

The local encoding equation of edge e G C^; in terms of its predecessors is 



Ve{x) = Dvp{x)+D 

peP{e)nCE 



Vpr (x) 

peP(e)nP(C) e'eP(e)nCB p'eP'{e)nP{C) 

(7) 

where P'(e) = {p e E \ end{p) = start{e)]. Note that P(e) C P'(e) but the 
converse is not true in general. For instance, in Example 4 P'{ei^) = {65,68,617} while 
P(6i3) = {65,617} (see Figure 6). 

The second sum in the formula brings the updated versions of the symbols that enter 
the knot at that point, while the double sum in the third term of the formula takes care 
of removing the old symbols. 

This results in the following local encoding equation of each edge in the cycle in terms 
of the predecessors of the cycle: 

v,{x)= Y D'^^''^r{p,e)Vp{x) (8) 
peP(c) 

One can see that the previous case is just a particular case of this one, since the line 
graph that will be associated to a simple flow cycle will always contain a simple cycle 
itself, and the corresponding transfer function between each p G P{C) and each edge e in 



the cycle will be r(p, e) = d{p, e) as was shown in 3.2.1 

A final observation at this point is that when the network is fiow acyclic or contains 
only simple fiow cycles, the global encoding equations will only contain polynomials on D, 
and not rational functions. Rational functions will be the result of using Mason's formula 
on knots. 

Considerations about how to decode will be discussed in the Appendix. 
We conclude the current Section by summarizing the complete LIFE* algorithm. 
ALGORITHM LIFE* 

Input: A directed multigraph G; a set of fiow paths /. 
Initializing: Vt : 

• Et = {e\e G start{e) = Si,i = 1, . . . , h}, 

• {ve{x) = Dai{x) = ai{x — l)\ e G Et, start{e) = Si} or, equivalently, Mf = DI^,, 
Vt G T. 

Main loop: Select an edge e for which the encoding equations have not yet been deter- 
mined, but for which the global encoding equations of all the predecessor edges in 
P{e) have been determined. Then proceed with the update of the set of current 



edges and current encoding equations as described in Subsection 3.1 
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If selecting such and edge is not possible, then a flow cycle has been encountered. 
Follow the procedure explained in Subsection 3.2 



Output: For each edge, the local encoding as given by (|2]) is produced. At the end of 
the algorithm, Et = {e \ end{e) = t} and Vt is still a set of h linearly independent 
equations from which t can recover the input. 

Examples of application of the algorithm can be found in the Appendix. 

4 Practical considerations 

The algorithm will execute the main loop at most \E\ times. The exact complexity 
depends on details of the algorithm not discussed here. However the complexity of the 
LIF and LIFE algorithms are similar. For discussions on the complexity of the algorithm 
we also refer the reader to [3]. 

The encoding presented here follows a flow path graph given for a network. This flow 
path graph is not necessarily unique and the choice made when computing the flow path 
graph determines much of the possible encodings that can be achieved. Which flow path 
graph is the best choice remains an open problem. First one should consider in which way 
the solution wants to be optimal (minimal delay, minimal number of link used, minimal 
number of encoding nodes ...). Some notions of minimality in the flow path graph can be 
considered that we will not discuss here. Also we will not discuss the different strategies 
that can be used in order to compute a flow path graph. 

Once the flow path graph for the network has been computed, the algorithm proceeds 
by following a topological order of the edges whenever that is possible (until a flow cycle 
or knot is found). However, this topological order is not unique. In many cases there is 
a certain choice to be made at each step about which edge will be encoded next of the 
several that follow in the order. This choice might in certain cases influence the total 
amount of delay necessary for the encoding. Examples can be given in which different 
choices of order lead to different final amounts of delay. Which ordering is most convenient 
for each flow path graph is also an open problem. 

Another consideration to take into account is that the presence of added delay means 
that the nodes at which the delay has to be introduced must have memory elements to 
store the symbols that have to be 'artificially' delayed. In most cases the sinks will need 
to use memory in order to be able to solve the equations. In any case the maximum delay 
used is finite. In case no extra delay needs to be added to the maximum delay needed on 
each path will correspond to the total length of that path from source to sink. 

The OSR and PI principles are also not necessary, but they help to keep the encoding 
simpler. Still, encodings can be found for flow cyclic networks in which the principles are 
not respected. Not respecting the PI principle is in fact equivalent to choosing a different 
flow graph path. 

4.1 Network preceding 

The inverse matrix of the encoding equation system may contain rational functions with 
denominators not on the form of D\ for some constant i. If so, the encoding is 'catas- 
trophic' in the sense that an error occurring in one of the transmissions can result in an 
infinite sequence of errors at the decoding sink. In order to avoid that, once the encoding 
has been computed using the LIFE* algorithm, one can compute the polynomial which is 
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maximum common divisor of all the rational functions resulting in the encoding process 
and introduce a pre-coding of the symbols generated by the sources, multiplying them by 
that maximum common divisor before they are introduced in the network. Alternatively, 
we can carry out this precoding locally in the nodes where a path enters a knot. We omit 
the details. 

After this precoding is introduced, the network code as viewed from the perspective 
of the sink is polynomial, and any error that might occur will cause only a limited error 
propagation that can be handled by a suitable error correcting or erasure restoring code. 

5 Conclusions 

The LIFE* algorithm is able to encode any given network with polynomial complexity 
and over the binary field. The addition of delay at some nodes is not a major drawback. 
In fact, any network encoding will in practice have intrinsic delay associated with it, and 
the delay will differ over the various paths. Thus in most cases, LIFE* does not need to 
introduce extra delay. In the few cases in which we actually need to introduce extra delay, 
this extra delay is what allows us to get the encoding on the binary field, which would 
have been impossible otherwise. For networks where the LIF/LIFE algorithms work, 
LIFE* will perform with essentially the same complexity as the others, i. e., there is no 
known more efficient algorithm in these cases. If the network contains knots, which many 
practical networks do, no other known algorithm works, but for LIFE* the complexity 
may become dominated by the calculation of Mason's formula. The greedy approach to 
finding the coding coefficients for each edge performs essentially as in the acyclic cases, 
also for knots. 

Appendix: Examples 

We will show here how the LIFE* algorithm will find encodings for the different types of 
networks shown in Figure 2. 

Example 1 

The (2,4) combination network is presented in Figure 2 a). It is known ([3],[H], [5]) to be 
a network which when extra delay is not used requires a finite field larger than F2. We 
show here how the LIFE* algorithm will work on the flow graph given in Figure 3. 

In order to better follow the progress of the algorithm, we have assigned labels 
Ci, . . . , ei8 to the edges in the network following a topological order. Also, for simplicity 
in the notation we have called the two sources A and B, and a{x) and b{x) are the binary 
symbols released by the sources at time x. 

The flow paths are represented with a code of colors and patterns in order to make 
is visually easy to follow. Each sink has a color assigned and each of the 2 sources has a 
pattern assigned (solid for A, dashed for B). The flow path from a source to a sink will 
be drawn in the color of the sink and with the pattern of the source. 

The LIFE* algorithm starts by setting the following initial values: 

Et,=Et, = --- = Et^ = {ei, 62} 
Vei{x) = Da{x) = a{x — 1) 
Ve2{x) = Db{x) = b{x — 1) 
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Figure 3: A flow path graph for the (2,4) combination network. 



D 



Hence 

\ D 

Now the algorithm enters the main loop: 



z = 1,. . . ,6. 



For encoding 63 we can observe that the only predecessor is ei, and thus 

Edge 63 is in the flows to sinks ti, t2 and ^3, so we update the corresponding sets of 
edges and matrices 

^ti = = = {63, 62} 



^i,3(ei)+2 Q 

D 



Clearly the choice ^63(61) = makes all the matrices non singular. Thus the encoding 
of 63 is 

Ve^{x) = Dve^{x) = D'^ a{x) = a{x — 2) 
64 has two predecessors, ei and 62- 

Ve^{x) = D'^^^^'^+he.ix) + D'^^^^^^+he,{x) = 0'^^'^^^^+^ a{x) + D'^^^^^^^%{x) 

Edge 64 is in the flows to sinks ti, and ts, so we update the corresponding sets of 
edges and matrices 

Et^ = {e3,e4},Et^ = {e4,e2},Et, = {64,63} 
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Again one can see that the choice ie^iei) — ieS^T) — makes all three matrices non 
singular. Hence 

Ve^{x) = Dve,{x) + Dve^ix) = D'^a{x) + D%{x) = a{x - 2) + b{x - 2) 

• Edge 65 has ei and 62 as predecessors, and the form of the encoding is 

Ve,{x) = D'^5^^'^+he,{x) + D'^^^^^^+'^Ve,{x) = D'^5^^'^+^a{x) + D'^^^^^^+^b{x) 

Edge 65 takes part in the flows to sinks ^2, ^4 and Iq, and the corresponding updating 
of edge sets and matrices is as follows: 

Et2 = {es, 65}, Et^ = {e4, 65}, Et^ = {e^, 62} 

Now clearly any value of ^65(61) and ^65(62) will make matrices M^^ and M^g non 
singular, but in order to get M^^ non singular we need those two values to be 
different, hence setting both equal to does not work in this case. A possible choice 
would be ^65(61) = 1, ^65(62) = 0, which gives us the next encoding. 

Ve^{x) = D^Ve^{x) + Dv^^^x) = D^a{x) + D%{x) = a{x - 3) + b{x - 2) 

• In the same manner we work with edge Cq, which has only one predecessor, namely 
62 and following the same procedure as before we can see that setting the only 
unknown exponent to will give a correct encoding. 

Veeix) = Dve^{x) = D%{x) = b{x - 2) 

Remcirk: 

The particular case we have seen in the encodings of edges 63 and ee, that is to say, an 
edge with only one predecessor is always solved in the same manner, copying the symbol 
carried by the predecessor and adding the natural delay unit. This corresponds to 

Ve{x) = DVp{x) 

when p is the only predecessor of e, that is to say,P(e) = {p} . (This means that the 
exponent ie{p) has been chosen to be 0.) 

The updated matrices will keep full rank since, for each t e T{e) the corresponding 
updated matrix will be the result of multiplying by D the elements of one of the columns 
of the old matrix, which does not alter the rank of the matrix. □ 
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The rest of the encoding steps in this example are trivial in that sense, since all the 
rest of the edges have an only predecessor. 

To complete the example we will illustrate how sinks can decode, this will also show 
what the delay means at the receiver end. 

Let us focus on sink tg. According to the encoding just computed this sink will receive 
at time x the symbols a{x — 4) + b{x — 3) and b{x — 3). 

Since we are assuming a{x) — b{x) — for all negative x, sink will receive zeros 
on both channels until time a; = 3 in which it receives + 6(0) and 6(0). This obviously 
gives him the knowledge only of symbol b{0). But at time a; = 4 it receives a(0) + 6(1) 
and 6(1). The knowledge of 6(1) allows it to recover a(0), which completes the recovery of 
the symbols of generation 0. Proceeding in the same way it will complete the recovering 
of the symbols of generation x at time x + A. The total delay observed by sink te is 4, 
which in this case coincides with the maximum power of D used in the encoding equations 
arriving in tg- 

This can also be interpreted in terms of matrices. 

f \ 

If we create a vector with the symbols that arrive at at time x and denote it as 
[?^t6,i(^)) ^t6,2(^)]) the encoding process can be described as 

[a{x),b{x)]Mt^ = [rt6,i(x),rt6,2(a;)] 

(which is equivalent to saying that rtf^,i{x) = a{x — 4) + b{x — 3),rt^^2{x) = b{x — 3)). 
Now the decoding process can be described as 

[a{x),b{x)] = K,i{^),rt,,2{x)]M-^'^ = K,i{x),rt,,2{x)]^ 

that is to say 

a{x) = rte,i{x + 4) + rte,2{x + 4) 
b{x) = rtg,2(a; + 3) 

which again shows how, to recover the symbols in generation x, sink has to wait 
until receiving symbols at time x + 4. 
Remcirk: 

In general the delay experienced by each sink is lower bounded by the maximum length 
of the flow paths arriving at it from the h different sources and upper bounded by the 
maximum power of the delay operator D used in the global encoding equations of the 
edges arriving at that sink. □ 

The upper bound is not always tight. To illustrate this let us consider the decoding 
that sink t4 in the example has to do. Despite the maximum power of D for that sink is 
4, it is easy to see that ^4 will complete the recovery of generation x at time x + 3. (But 
in addition it is absolutely necessary for ^4 to keep one memory element in order to be 
able to decode.) 



IS 
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Example 2 




Figure 4: The unique flow path graph for Example 2. 

Example 2 

This example shows a network which is link cyclic but flow acyclic. The encoding process 
will work analogous to what was shown in the previous example. 

The edges in Figure 4 have been numbered according to a topological order and one 
can observe that each edge has in fact only one predecessor, hence the encoding becomes 
trivial. We simply show here the result obtained. 

Ve^ix) = a{x — 1) 

Ve^{x) = b{x - 1) 

a{x — 2) 
b{x - 2) 
a{x — 3) 
b{x - 3) 

As we can see, link cyclic but flow acyclic networks do not present any additional 
problem for encoding, they behave exactly as the acyclic networks did. 

Example 3 

Here we deal with a flow cyclic network that contains a simple flow cycle. Figure 5 shows 
the unique flow path graph for this network. 

The initialization will give us the encoding of the flrst 9 edges 

Ve-,{x) = Ve^ix) = f eg (x) = a{x ~ 1) 
Ve^ix) = Ver^{x) = f eg (x) = 6(x - 1) 
Verix) = Veg{x) = f eg (x) = c{x - 1) 

Now no more edges can be visited following a topological order. The cycle has set of 
edges Ce = {610,611,612} and set of predecessors -P(C) = {61,65,69}. 
The structure of the local encoding equation in the cycle will be 

vc{x) = D^^(^iV(6i,C)t;e,(x) + D^^(^^)r(65,C)t;e5(a;) + /^^^(^«V(69,C)t;eg(x) 

Since we are in the case of a simple flow cycle, we can follow the formula given in |3.2.1 
for the local encoding equation of each of the three edges in Eq- 



VeAx) = 
Ve,{x) = 
Ve.ix) = 
Veeix) = 
Verix) = 
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Example 3 




Figure 5: The unique flow path graph for Example 3. 



Ve^{x) = Z}*^(^i)+^("i'"^\ei(a;) + D'^^'''^+'^^^^^''^^Ve,{x) + D'^^^'^+'^^^'''^^Ve,{x), J = 10, 11, 12 

Inspection of the graph shows that 

d{ei,eio) = l,rf(e5,eio) = 3,c?(e9,eio) = 2 
d{ei, en) = 2, d^e^, en) = 1, d{eg, en) = 3 
d{ei, 612) = 3, (i(e5, e^) = 2, d{eg, 612) = 1 

and the substitution in the above expression gives 

feio(x) = D'ci^^^+^a{x - 1) + D'c^^^'''^+^{x - 1) + D'cM+'^c{x - 1) 
VeSix) = D'c(^^)+^a{x - 1) + D'c<^''^+^b{x - 1) + D'c<^^'^-^^c{x - 1) 

VeSix) = D'c(^^)+'^a{x - 1) + D'C<^^'^+%{X - 1) + D*C-(e9) + lc(a; _ 1) 

The full rank invariant condition must be checked for edge eio in the flow to sink ^3, 
for edge en in the flow to sink ti and for edge ei2 in the flow to sink t2. This gives us the 
following matrices: 

/ £,ic(ei)+3 \ / D L)^c(ei)+4 Q \ /DO L)^c(ei)+2 \ 

Mi, = L)*c{e5)+2 ^ , = D*C-(e5)+3 Q , = D L)^c(e5)+4 

\ D / \ D*C?(^9)+2 jj J \ L)ic(e9)+3 y 

Clearly any value of «c(ei), ^0(65) and ici^g) satisfles the full rank invariant and we 
choose the simplest one setting the three unknowns to be 0. 

The global encoding equations of the edges in the cycle are as follows 





[x) 


= a{x 


- 2) + b{x 


-4)- 


- c(x 


-3) 




[x) 


= a{x 


- 3) + blx 


-2)- 


- c(x 


-4) 


^ei2 


[x) 


= a{x 


- 4) + b{x 


-3)- 


- c{x 


-2) 
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The encoding now of edges 613, 614 and 615 is trivial since each of them has only one 
predecessor. 

Veiaix) — a{x — 4) + b{x — 3) + c{x — 5) 
Vei^{x) — a{x — 5) + b{x — 4) + c{x — 3) 
Veig{x) — a{x — 3) + b{x — 5) + c{x — 4) 

The delay at the final receivers is 4, even when the maximum exponent of D in the 
equations received by the sinks is 5. Besides there is some extra memory needed in order 
to decode. For instance, recovering the element a(l) and hence completing the generation 
1, can be done by sink ti at time x = 5, provided it kept in memory the element c(0). 

A shght modification could be done for the encoding of the edges whose predecessors 
lie in the cycle, in such a way that they get the last updated values, for instance, edge 613 
can benefit from the fact that edge eg enters in the same node from which 613 exits and 
hence get an updated version of the symbol carried by eg, then the encoding of ei3 would 
be 

Veisix) — a{x — 4) + b{x — 3) + c{x — 2) 

which is actually the same encoding that has the edge ei2, and results in smaller memory 
needed at the receiver ti. 
Remcirk: 

In general, using this last observation, the local encoding equation of an edge e ^ Ce 
with a predecessor pc G -P'(e) H P{C) would be 

Ve{x) = D'-^'^^Vsc{x) + ^ D'^^^^Dvp{x) 

peP(e)\{pc} 

where sc is the successor of pc that lies in the cycle, that is to say, the edge in 
with start{sc) — start (e) = end{pc)- 

We finally remark that in a simple flow cycle, the element sc is are unique, even in 
there are several elements pc in P'{e) n P{C). □ 



Example 4 

In this example we show how to work with a knot. Figure 6 shows the essentially unique 
flow path graph for the network given. 
The initialization values are 

VeAx) = Ve2ix) = Ve.^{x) = a{x - 1) 
Ve^ix) = Ve.ix) = Vee{x) = b{x - 1) 
Ve.,{x) = Ves{x) = Ve^ix) = c{x - 1) 
Ve-,o{x) = 'yeii(x) = Ve:,^{x) = d{x - 1) 

No more edges can be visited in topological order because Ce — {ei3, eu, ei5, Ciq, en} 
form a flow cycle. In fact it is a non simple cycle, since it contains two fiow cycles 
ei3 -< ei5 -< ei7 -< 613 and 614 -< eig -< en -< eu, both sharing the edge eiy. Hence we are 
in presence of a flow knot. (See Figure 7 a)) 

The predecessors of the knot are Pc = {62, 65, eg, en}. 

The general structure of the local encoding equation for the knot is 

vc{x) = D'^^''^T{e2, C)ve,{x)+D'^^'''^T{e5, C)ve,{x)+D'^^"'^T{es, C)ve,{x)+D''^^''''^T{eu, C) 
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gure 6: A flow path graph for the network in Example 4. 
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In order to find the local encoding of each edge in the knot we have to use Mason's 
formula on a line graph to compute the transfer function for each of the symbols carried 
by the predecessors of the knot. 

The line graph is shown in Figure 7 b). The nodes correspond to the edges in Cg. 
The edge connecting eis to eis is drawn because of the flow path from source B to sink 
^2, the edge connecting eis to 617 is determined by the flow from source A to sink ti. In 
the same way we draw the other connections in the line graph following the flow paths. 

The branch gain of each connection is D. Mason's formula is as follows 

r{ej,ek) = '—^ 

where A = l + ^Cj — ^ CjCj + ■ ■ ■, Fi{ej, e^) is the function corresponding to the i — th 
forward path form Cj to and Aj(ej, e^) is defined as A but counting only the cycles in 
the circuit that are disjoint with the i-th forward path. Here we are using the notations 
in and we refer the reader there for a more detailed explanation of Mason's formula. 

• We will now focus on the symbol that enters through edge 62 into 615. Its itinerary 
through the knot is shown in Figure 8 (a). 





Figure 8: The itinerary of symbols carried by 62 in the knot of Example 4. 

The edges used by the flow path that carries that symbol are represented as bold 
lines, but we observe also that when the symbol arrives at node end{en) this node 
will distribute it not only to edge 614, but since there is a flow path connecting edge 
ei7 with edge 613, the symbol in question will travel also on edge 613 and in the 
same way we can see that it will travel also on edge Cig. That is represented in 
dashed lines in Figure 7 (a). On the other hand, using memory, the node end{e2) = 
end{ei3) = start{ei^) can remove the contribution of the old symbol that arrives 
back at it through 613, so we can then remove the connection between edges 613 and 
ei5 from the line graph. Hence the actual line graph followed by the symbols that 
enter the knot through edge 62 is shown in Figure 8 (b). 
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In that graph there is only one cycle, which has length 3, namely {en, eu, eie, 617}. 
All the transfer functions will have as denominator the function A — 1 + D^. 

For the transfer function corresponding to edge 615, Figure 9 (a) shows that the only 
forward path is node disjoint with the only cycle of the graph, hence Fi{e2, 615) = D 

and Ai(ei5, 615) = 1 + Finally r(e2, 615) = ^^^^^ ^ = D which is the expectable 
result. 




Figure 9: Computation of Mason's formula from edge 62- 

Now we compute the transfer function corresponding to edge eiy. Figure 9 (b) shows 
the only forward path, which is not node disjoint with the only cycle of the graph, 
hence ^1(62,617) = D"^ and Ai(e2,ei7) = 1. Finally r{e2,en) = j^^. 

In the same way all the other transfer functions can be computed (drawings of the 
corresponding forward paths can be seen in Figures 9 c) to e)). 

T(e2,ei4) = 7— — =^,r(e2,ei6) = . , „^ ,r(e2,ei3) = 



1 + ^3' V 1 + D3' V 1 + ^3 

Specially interesting is the transfer function corresponding to edge 613 since it gives 
the function of the old symbol that has to be removed when passing again through 
node end{e2) — end{eis) — start{ei^). 

• If we now focus on the circulation in the knot of the symbol carried by edge 65 we 

observe (Figure 10) that the flow path will transport it through edges 613 and 615 
(bold line in the flgurc) and node en(i(ei5) will send it back to node end{e^) through 
edge ei7 (dashed in the figure), but node end{e^) will not send that symbol on edge 
ei4 and hence the symbol does not travel along the two cycles in the knot, but only 
on the cycle {ei3, 615, 617, 613}. The computation of the transfer functions is then 
straightforward . 

T(e5, 613) = D, r(e5, 615) = r(e5, 617) = D^, T(e5, eu) = 0, r(e5, eie) = 
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Figure 10: Itinerary of symbols carried by 65 in the knot of Example 4. 



• The symbol entering the knot via edge eg follows a similar trajectory to that of the 
one entering via edge 65. 

r{es: eu) = -D, T{es, eie) = -D^ T(e8, 617) = D^, T{es, 613) = 0, T{es, 615) = 

• Finally, the symbol entering via edge en follows an itinerary identical (considering 
symmetry) to that entering via 62 already studied. 

/ \ D-(l+D^) 7-, / X D^-l 

T(eii, eie) = ^^j^s = -D, T(eii, 617) = 

''■(,eii,ei3j = i:^,r[eii,ei5) = jqios > '''leii, euj = Y^r^s 

Now that all the transfer functions have bee computed one can use formula in Equation 
(5) to compute the encoding of each edge in the knot. 

We will go in detail with the computation of the encoding of edge 613. 

Ve,,{x) = D^^(^^h{e2, ei3)v2(a;) + D^^^'^h{e5, ei3)v5ix) + 

= D'oi^^^j^a{x - 1) + D'^^^^Wb{x - 1) + D'c^^^^Oc{x - 1) + D^^^^") - 

In the same way the encoding of the other edge in the knot can be computed using 
Equation (6). 

Ve,^{x) = D'^'^^-'^j^aix - 1) + L)*c'(e5)o5(a; _ 1) + D^'^^^^'^ Dc{x - 1) + D'^^^^^'^ ^d{x - 1 
Ve^^{x) = D'c^^^Wa{x - 1) + D'(^^^^W%{x - 1) + D'^^^^^Oc{x - 1) + D'^i^^^^ j^d{x - 1) 
= D'^^^^^j^a{x - 1) + D'(^^^^^Ob{x - 1) + D'<^^^^W^c{x - 1) + D'^^^^^^Dd{x - 1) 

Ve^,{x) = D'c(^^^j^a{x - 1) + D'c(^5)D%{x - 1) + D'ci^^^D^c{x - 1) + D'ci^^^^^d{x - 

Now we will show how each edge can compute its encoding based on its predecessors 
and using formula in Equation (5). Again we will go in detail with edge 613. 
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P{ei3) = {65, en}, P{en) n P{C) = {eg}, P(ei3) n = {eir} 
/"(ei3) = {es, eg, e^}, /"(eia) n P(C) = {es, eg} 

The direct application of formula in Equation (5) to this case gives the following 

= D (^D'c(e2)^a(a; - 1) + D'^^^^^D%{x - 1) 

+ D'ci^^W^c{x - 1) + D'^'^^^^^ j^d{x - 1) 
+D'c^''^^Db{x - 1) 

-D^^^^^WD%{x - 1) - D^c^^^WD^cix - 1) 



= L)ic(e2)_^a(^ - 1) + D'c^^^^Dh{x - 1) + L>^'^(^«)Oc(x - 1) + ^^^^(a; - 

In a similar way all the other edges can get their encoding using those of its predecessors 
and formula in Equation (5). 

Next point is determining the unknowns ic(e) for each edge in C by checking the 
full rank conditions. To be precise, the full rank condition must be checked for edge 
ei3 in the flow to sink ti, for edge ei4 in the flow to sink to^ and for edges 615 and Cig 
in the flow to sink t2- A careful exam of the corresponding matrices will show that 
ic{&2) = ic{&b) = ic{es) = ic(eii) = is a valid choice. 

Finally, a similar remark to the one made in Example 3 gives us Vei^{x) — Ve^^{x) and 
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Also we have Ve-^,^{x) 



DVe,,{x),Ve,„{x) 



The flnal encoding matrices at the sinks are 
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As can be seen the elements in the matrices are now rational functions (typical case af- 
ter traversing a knot) . This however represents no extra difficulty at the sink. For instance, 

if we focus on sink t2, the received symbols at that sink are Tt2,i{x) — Vi{x),rt2,2{x) — 
Vi(j{x)^rt2,?,{x) = V2o{x) and ^^2,4(2;) = ^12(2;). Their expressions in terms of the source 
symbols are given in each column of the matrix Mf^- In order to retrieve the source 
symbols t2 will multiply the received ones by the inverse of matrix M^^. 



[a{x),b{x), c{x), d{x)] 



[rt,,i{x),rt,,2{x),rt,,3{x),rt,,4{x)]Mt^^ 
[rt2,iix + l),rt2^iix + 2) + rt.^^2ix + A) + j^rt^,4ix), 



T^rt2^i{x) +ri2,3(a; + 4) + rt2,4ix + 2),rt^^4{x + 1)] 

The way to deal with expressions like ^^^3 ^t2,i(^) ^1 to keep in memory three 
local variables that we will call rt2,i,i{x). They will all be initialized as 0, and at time x 
one of them will be updated and the other two keep the same as follows: 

n2,i,iAx) = n2,i,i.(a^-i) + ^*2,i(a^) 

rt2,i,ii^) = n2,i,ii^ - 1) for ^ e {0> 1> 2} \ {i^} 

where x — Sqx + ix with G {0, 1, 2}, that is to say, ix is the remainder of the integer 
division of x by 3. 

In this way -^^^3 i"t2,i{^) = ^t2,i,i(a;-i) ~ 1)' receiver t2 does not need to keep 

an infinite memory, despite the aspect of the equations received. 
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Example 5 



We will briefly show here one more example in which a more complicated network is dealt 
with. 




Figure 11: The network of Example 5. 



The network is shown in Figure 11. It has six unit rate sources, labeled A, B, . . . , F. 
and three sinks, labeled ti,t2,^3- 

The other nodes have been labeled 1,. . . , 7. 




Figure 12: Flow paths in the network of Example 5. 

The flow path graph for such a network is (essentially) unique. The flow paths starting 
in sources A and B are shown in Figure 12 in solid bold and dashed bold lines respectively. 
The flow paths from the other sources are the same but with a rotation of 120 degrees to 
the right or to the left. 

This implies that the flow path graph of that network contains a knot formed by 6 
simple cycles, each of length 3. Figure 13 a) shows the knot while Figure 13 b) shows the 
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(a) (b) 



Figure 13: The knot and its line graph in the network of Example 5. 

corresponding line graph constructed following the flow paths through the knot. 

We will focus now on the way the symbols are carried by the edge that connects source 
A with node 1 travel through the node. Mason's formula will be used to compute the 
transfer functions from that edge (which is one of the predecessors of the knot) to any 
other edge in the knot. Let us call a that edge, that is, a — {A, 1) e P{C). In Figure 14 

a) we show the itinerary of the symbol from a. The bold solid lines are the actual flow 
path and the dashed lines are the edges not belonging to the flow path of the symbol 
but that will nevertheless carry instances of that symbol due to the connections in the 
knot. We can see that the only edge in which that symbol do not travel is en. Figure 14 

b) shows the corresponding modified line graph that allow us to compute the transfer 
function t(q;, 12). The bold dashed arrows show the position of the 4 simple cycles in 
that graph and the two bold lines (one black and the other grey) show the trajectories 
of the two different forward paths from a to ei2. Direct application of Mason's formula 
gives 

~ 1 + 4L>3 + SL'e ~ 1+^03 
The rest of the transfer functions from a to edges in C are 

t(q;, ei) = D, r(a,e2) = i:*^ + r(a,e3) = j^, r(a,e4) = Y^, 

T(a,e5) = Y^, r(a,e6) = Y^, ^(",67) = ^^, T{a,es) = 

'r(a,e9) = j^, T(a,eio) = Y^, r(a,eii) = 0. 

The transfer functions from the other predecessors of the knot to the edges in the knot 
arc analogous. In fact, they can be derived from the ones already computed by simply 
taking into account the multiple symmetries that this network presents. 
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Figure 14: Itinerary of the symbol from edge a through the knot in the network of 
Example 5. 



If we denote by /3 the edge that connects source B with node 1, 7 the one that connects 
C with 5, 5 the one that connects D with 5, and finally e and the edges connecting 
sources E and F respectively to node 3, following the procedure of LIFE* we obtain the 
following global encoding equation for edge 65 



(/3). 



n8 r)8 



1 + V ^ ' 1 + V ^ ' 1 + 

Encoding for other edges in the network will be analogous. 

The full rank condition will be satisfied when choosing ec{cy) — e.c{(3) — 6^(7) = 
ec(5) = ec(e) = ec(0) = 0. 

The local encoding equation of edge 65 in terms of its predecessors is as follows. 

Ve^{x) = Dve^{x) + DD'^^^^v,{x)) - D (L>'c(^V(e, 64)^,(0;) + 64)^^(0;)) 
The decoding matrix for sink ti has the form 



d{x)+D'^^'^D'^e{x) 



( 



M- 



1 



D 



1 

2 , 



1 

1+D6 

D2 



^ + 1+D6 
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